Dynamic scaling and aging phenomena in short-range Ising spin glass: 
Cuq 5Coo.5Cl2-FeCl3 graphite bi-intercalation compound 
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Static and dynamic behavior of short-range Ising-spin glass Cuo.5Coo.5Cl2-FeCl3 graphite bi- 
intercalation compounds (GBIC) has been studied with SQUID DC and AC magnetic susceptibility. 
The T dependence of the zero-field relaxation time r above a spin-freezing temperature Tg (= 3.92 
± 0.11 K) is well described by critical slowing down. The absorption x" below Tg decreases with 
increasing angular frequency uj, which is in contrast to the case of 3D Ising spin glass. The dynamic 
freezing temperature Tf{H,uj) at which dMFc{T, H)/dH — x{T,H — 0,tj), is determined as a 
function of frequency (0.01 Hz < cj/27r < 1 kHz) and magnetic field (0 < -ff < 5 kOe). The dynamic 
scaling analysis of the relaxation time t{T,H) defined as r = l/uj at T = Tf{H,Lj) suggests the 
absence of SG phase in the presence of H (at least above 100 Oe). Dynamic scaling analysis of 
x"{T,Lu) and r{T,H) near Tg leads to the critical exponents (/3 = 0.36 ± 0.03, 7 = 3.5 ± 0.4, u = 
1.4 ± 0.2, z = 6.6 ± 1.2, ^ = 0.24 ± 0.02, and 9 = 0.13 ± 0.02). The aging phenomenon is studied 
through the absorption x"{^ji) below Tg. It obeys a (tjt)"'' power-law decay with an exponent 
b" ~ 0.15 — 0.2. The rejuvenation effect is also observed under sufficiently large (temperature and 
magnetic-field) perturbations. 

PACS numbers: 75.50.Lk, 75.40.Gb, 75.40.Cx 



I. INTRODUCTION 

In recent years there has been much effort to under- 
stand phase transitions of three-dimensional (3D) Ising 
spin glasses (SG's) with short-range interactions. The 
dynamic behavior of the low temperature SG phase has 
been a subject of some controversy. There are basically 
two different pictures of the SG phase, the mean-field pic- 
ture (including Monte Carlo simulations)^ and the scal- 
ing picturei^i^i^i^i^''' 

The mean-field picture is originated from the replica- 
symmetry-breaking (RSB) solution of Parisi for the 
infinite-range Sherrington-Kirkpatrick (SK) modeli In 
this picture the SG phase has a complicated free energy 
structure with a number of degenerate minimum sepa- 
rated with each other by large free energy barriers. The 
SG phase survives in the presence of a magnetic field. A 
de Almeida- Thouless (AT) line separates the SG phase 
from a paramagnetic phase.® 

The scaling picture (so-called the droplet model) is 
based on renormalization group arguments In 
this picture there are only two thermodynamic states re- 
lated to each other by a global spin flip, and the im- 
portant excitations at low temperatures are droplets of 
overturned spins. In three dimensions a true phase tran- 
sition should exist only in zero field and no irreversibility 
should be seen in the presence of any finite magnetic field. 

Because of extremely long relaxation time, there is 
no direct measurement available that can demonstrate 
whether the SG phase exists even in zero magnetic 
ggj(^i9ji0jii rpj^g dynamic scaling analysis is required to 
overcome such an experimental difficulty in confirming 
the existence of thermally equilibrium SG phase in a fi- 
nite magnetic field. This problem has been addressed by 
Mattsson et alt? They have studied the dynamic suscep- 



tibility of a 3D Ising SG, Feo.sMno.sTiOa. The in-field 
scaling behavior of the dynamic susceptibility, which is 
remarkably different from the zero-field behavior, is ex- 
plained as demonstrating the absence of SG phase in a 
finite field. 

In this paper we study the slow dynamics of the 
SG phase in a short-range 3D Ising SG, Cuo.5Coo.5Gl2- 
FeCla graphite bi-intercalation compound (GBIC), using 
SQUID DC and AC magnetic susceptibility. This com- 
pound undergoes a SG transition at a spin freezing tem- 
perature Tg (= 3.92 ± 0.11 K) at zero magnetic field. 
Following a method used by Mattsson et al.;^ we de- 
termine the dynamic freezing temperature Tf(H,io) at 
which dMFc{T,H)/dH = x'{T,H = Q,uj), as a func- 
tion of frequency (0.01 Hz < w/27r < 1 kHz) and field (0 
< H <5 kOe) . The dynamic scaling analysis of the relax- 
ation time r(r, H) defined as r = 1/w at T = Tf{H,uj), 
suggests that the SG phase does not survive in the pres- 
ence of H (at least above 100 Oe). The critical expo- 
nents are determined from the dynamic scaling analy- 
sis for x"(r, cj) and t{T,H) and compared to those for 
Feo.sMno.sTiOa. The aging phenomenon is also studied 
through the time dependence of the absorption x"(a;,i) 
below Tg. It obeys a {ujt)~^ power- law decay with an 
exponent b" w 0.15 - 0.2. 

Cuo.5Coo.5Cl2-FeCl3 GBIC has a unique layered struc- 
ture where the Cuo.5Coo.5Cl2 intercalate layer (= Ii) 
and FeGla intercalate layers {—I2) alternate with a sin- 
gle graphite layer (G), forming a stacking sequence (- 
G-h-G-h-G-h-G-h-G-) along the c axis.i^-ia in the 
Cuo.5Coo.5Cl2 intercalate layer two kinds of magnetic 
ions (Cu^+ and Go^"*") are randomly distributed on the 
triangular lattice. The spin order in the Guo.5Coo.5Cl2 
layers is coupled with that in the FeCls intercalate layer 
through an interplanar exchange interaction, leading to 
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the spin frustration effect. 

The XY and Ising character of the present system are 
due to the easy-plane type anisotropy field iJ™* and the 
in-plane anisotropy field H™, respectively: 77^"* 3> -ff™- 
Although the in-plane structure of Cuo.5Coo.5Cl2 layers 
is incommensurate with that of graphene sheets, Cu and 
Co atoms tend to sit over hexagon centers of graphene 
sheets because of the interplanar structural correlation. 
The sixfold symmetry of graphen sheets imposes the in- 
plane spin anisotropy in the c plane. Because of this, 
the spins tend to align along the easy axis with sixfold 
symmetry in the c planed 

It is well known that the intercalate layers are formed 
of small islands in acceptor graphite intercalation com- 
pounds (GIC's).^^ In stage-2 C0CI2 GIC, for example, 
there is a charge transfer from graphene sheets to the 
C0CI2 intercalate layer. The periphery of such islands 
provides sites for charges transferred. The size of islands 
is on the order of 400 A. In spite of few structural studies, 
it is predicted that similar island structures exist both in 
the Cuo.5Coo.5Cl2 and FeCla layers of our system. 



II. EXPERIMENTAL PROCEDURE 

A sample of stage-2 Cuo.5Coo.5Cl2 graphite interca- 
lation compound (GIC) as a starting material was pre- 
pared from single crystal kish graphite (SCKG) by va- 
por reaction of anhydrated Cuo.5Coo.5Cl2 in a chlo- 
rine atmosphere with a gas pressure of « 740 Torr>i^ii& 
The reaction was continued at 500 °C for three weeks. 
The sample of Cuo.5Coo.5Cl2-FeCl3 GBIC was pre- 
pared by a sequential intercalation method: the in- 
tercalant FeCla was intercalated into empty graphite 
galleries of stage-2 Cuo.5Coo.5Cl2 GIC. A mixture of 
well-defined stage-2 Cuo.5Coo.5Cl2 GIC based on SCKG 
and single-crystal FeCla was sealed in vacuum inside 
Pyrex glass tubing, and was kept at 330 °C for two 
weeks. The stoichiometry of the sample is represented 
by Cm(Cuo.5Coo.5Cl2)i-c(FeCl3)c. The concentration 
of C and Fe (m and c) was determined from weight 
uptake measurement and electron microprobe measure- 
ment [using a scanning electron microscope (Model Hi- 
tachi S-450)]: m = 5.26 ± 0.05 and c = 0.53 ± 0.03. 
The (OOi) x-ray diffraction measurements of stage-2 
Cuo.5Coo.5Cl2 GIC and Cuo.5Coo.5Cl2-FeCl3 GBIC were 
made at 300 K by using a Huber double circle diffrac- 
tometer with a MoKa x-ray radiation source (1.5 kW). 
The c-axis repeat distance of stage-2 Cuo.5Coo.5Cl2 GIC 
and Cuo.5Coo.5Cl2-FeCl3 GBIC was determined as 12.83 
± 0.05 A and 18.81 ± 0.05 1, respectively 

The DC magnetization and AC susceptibility were 
measured using a SQUID magnetometer (Quantum De- 
sign, MPMS XL-5) with an ultra low field capability op- 
tion. First a remnant magnetic field was reduced to zero 
field (exactly less than 3 mOe) at 298 K for both DC mag- 
netization and AC susceptibility measurements. Then 
the sample was cooled from 298 K to 1.9 K in a zero 



field, (i) Measurements of the zero-field cooled suscep- 
tibility (xzFc) md the field- cooled susceptibility (xfc)- 
After an external magnetic field H {Q < H <\ kOe) was 
applied along the c plane (basal plane of graphene layer) 
at 1.9 K, xzFC was measured with increasing temper- 
ature (T) from 1.9 to 20 K. After annealing of sample 
for 10 minutes at 50 K in the presence of iJ, xfc was 
measured with decreasing T from 20 to 1.9 K. (ii) AC 
susceptibility measurement. The frequency (/), magnetic 
field, and temperature dependence of the dispersion (x') 
and absorption (x") was measured between 1.9 to 20 K, 
where the frequency of the AC field is / = 0.01 - 1000 
Hz and the amplitude h is typically h = 50 mOe. 



III. RESULT 

A. xzFc{T,H) and XFc{T,H) 

Figures H^a) and (b) show the T dependence of 
the ZFC and FC susceptibilities {xzfc and xfc) of 
Cuo.5Coo.5Cl2-FeCl3 GBIC at various H, where the mag- 
netic field (H) is applied along the c plane perpendicular 
to the c axis. The T dependence of xzfc and xfc is sim- 
ilar to that of Feo.5Mno.5Ti03r^^ The susceptibility xzfc 
at H = 1 Oe exhibits a cusp around 4.5 K , while xfc is 
nearly temperature independent below 3.9 K. Note that 
Xfc along the c axis at H = 1 Oe is on the order of 
one-tenth of xfc along the c plane around 3.9 K. Figure 
Eta) shows the T dependence of the difference S defined 
by i5 = Xfc ~ Xzfc at various which is the measure 
for the irreversibility of the susceptibility. The deviation 
of Xzfc from xfc at H = 1 Oe starts to occur below 10 
K, and drastically increases with decreasing T below the 
cusp temperature of xzfc- The infiection-point temper- 
ature of (5 vs r corresponds to a local minimum temper- 
ature of dS/dT vs T shown in Fig.|2tb). The inflection- 
point temperature of 6 drastically decreases with increas- 
ing H, leading to the AT line in the {H, T) diagram (see 
Sec. IIII ni for further discussion). The non- vanishing 5 
well above the inflection-point temperature is partly due 
to the islandic nature of Cuo.5Coo.5Cl2 and FeCls lay- 
ers. The growth of the in-plane spin correlation length is 
partly limited by the existence of small islands, leading 
to the smearing of the inflection-point temperature over 
the system. 



B. Relaxation time: x'(T,uj) and x"(T,lj) 

In Fig. Ola) we show the T dependence of the disper- 
sion x'{T,uj) at various frequencies (/ = 0.01 — 1000 Hz), 
where lu {— 27r/) is the angular frequency and h {— 50 
mOe) is the amplitude of the AC field. In Fig. |2Ib), the 
T dependence of xzfc and xfc at = 1 Oe is com- 
pared to typical data of x'(T, uj). The susceptibility xfc 
may correspond to the equilibrium susceptibility (xeq)- 
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FIG. 1: (a) and (b) T dependence of Xzfc and X-FC at various 
H for Cuo.5Coo.5Cl2-FeCl3 GBIC. H is applied along the c 
plane (graphene basal plane) perpendicular to the c axis. 

In Fig. Efa) we show the T dependence of the absorp- 
tion x"(T, w) at various /. Similar behaviors in x'(T, a;) 
and x"iT,ijj) have been observed in Fe-C nanoparticles 
by Hansen et al,^- Figure^Jb) shows the T dependence of 
the derivative dx"{T,uj)/dT. If we assume conventional 
critical slowing down on approaching the SG transition 
temperature Tg from the high-T side, the T dependence 
of the relaxation time r can be described hy^ 

r = r*(T/T,-l)-^ (1) 

with X = vz, where z is the dynamic critical exponent, v 
is the critical exponent of the spin correlation length ^, 
and r* is the characteristic time. In the analysis of r vs 
T, T corresponds to the local minimum temperature of 
dx" /dT vs T, or the peak temperatures of x" vs T and 
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FIG. 2: T dependence of (a) 5 Xfc-Xzfc) and (b) dS/dT 
at various H. The value of 5 is derived from Fig.Q The solid 
lines are guides to the eyes. 

x' vs T, where r is set equal to w^^. In Fig. we show 
the T dependence of r thus obtained. The least squares 
fits of the data of r vs T yield x = 10.3 ± 0.7, Tg = 
3.92 ± 0.11 K, T* = (5.27 ±0.07) x 10"^ sec for the local 
minimum temperature of dx" /dT vs T, x = 15.5 ± 1.5, 
Tg = 2.98 ± 0.18 K, T* = (8.56 ± 0.05) x 10"^ sec for 
the peak temperature of x" vs T, and x — 15.3 ± 2.0, 
Tg = 3.57 ± 0.27 K, r* = (3.08 ± 0.10) x 10"^ sec for 
the peak temperature of x' vs T. The parameters thus 
obtained are rather different. The relaxation time r* is 
much larger than a microscopic relaxation time tq (typ- 
ically 10"^" sec). Here we assume that the local mini- 
mum temperature of dx" /dT vs T at a; = corresponds 
to the spin freezing temperature Tg. Tg is close to a tem- 
perature below which xfc at _ff = 1 Oe becomes nearly 
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FIG. 3: (a) T dependence of x' at various /. /i = 50 mOe. H piG. 4: T dependence of (a) x" and (b) dx"/dr at various 
= 0. (b) r dependence of xzfc and xfc {H = 1 Oe), which j /j = 5Q mOe. = 0. 
is compared with that of x' at / = 0.01, 1 and 1000 Hz (H 
= 0). 



constant. The value of x (= 10.3 ± 0.7) is rather close 
to that predicted by Ogielski for the 3D ±J Ising SG 
(x = 7.9 ± 1.0)2a from Monte Carlo (MC) simulations. 
The situation is a little different for Feo.sMno.sTiOaii^ 
In the analysis of t vs T, T is determined either as the 
maximum of x' or as the inflection point of x" ■ In both 
cases, T is well described by Eq. (P) with Tg = 20.95 ±0.1 
K, X = 10.0 ± 1.0, and logioTo = -12.8 ± 1.0, where tq 
is a microscopic relaxation time. The value of x is in 
excellent agreement with our result. 

Here we assume a generalized form of the relaxation 
time, more suitable to the description of SG behavior 
near T^iSi The relaxation time t(/,T), which is needed 
to flip the Z-sized cluster of spins, is governed by thermal 



activation over a barrier Bt(1), in such a way that 



T{l,T)^Tol'exp[BT{l)/kBT], (2) 



where I is in units of the lattice constant a. The energy 
barrier is described by Bt{1) = F(T)/'^, where ip is the 
barrier exponent. The wall stiffness Y{T) should vanish 
above Tg like in ferromagnets and is described by Y(T) — 
YQ{1 — T/Tg)^'^ below Tg, where v is the critical exponent 
of spin correlation length ^ and Yq is on the order of ksTg. 
Above Tg the exponential term in Eq.Q is nearly equal 
to 1, leading to the expression for the critical slowing 
down above Tg: t^{1,T) « tqI^. 
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FIG. 5: T dependence of relaxation time t, which is deter- 
mined from peak temperatures of x" vs T and x' "^s T, and 
local minimum temperature of dx" /AT vs T. The solid lines 
are fits to Eq. (0 (critical slowing down) with parameters x, 
Tg, and r* given in the text. 



C. Dynamic scaling of Tx"{T,iu) 

FiguresEfa) and (b) show the / dependence of x'{T, uj) 
and x"{T,uj) at various T in the vicinity of Tg (= 3.92 
K), respectively. The absorption x"(T, oj) curves ex- 
hibit different characteristics depending on T. Above Tg, 
x"{T, Lj) shows a peak at a characteristic frequency, shift- 
ing to the high /-side as T increases. Below Tg, x"{T,uj) 
shows no peak for / > 0.01 Hz. It decreases with in- 
creasing /, following a power law uj~°' . The exponent 
a" is almost independent of T: a" = 0.096 ± 0.003. Ac- 
cording to the fluctuation-dissipation theorem, the mag- 
netic fluctuation spectrum S{uj) is related to x"iT,Lu) by 
S{T,uj) = 2kBTx"{T,uj)/ujM Then S{T,lj) is propor- 
tional to \ which is similar to 1/lu character of 
typical SG. In contrast, x'C^^^) decreases with increas- 
ing / at any T near Tg: x' exhibits a power law . The 
exponent a' is weakly dependent on T: a' = 0.088±0.001 
at T = 3 K and a' = O.llliO.OOl at T = 3.8 K. The value 
of a' agrees well with that of a" . In fact, x" is related 
to x' through a so called "7r/2 rule": x" = -(7r/2)dx'/d 
Incj (Kramers-Kronig relation), leading to the relation 
a' = a" . Here we note that the observed / dependence of 
x" in Cuo.5Coo.5Cl2-FeCl3 GBIC is different from that in 
conventional spin glasses such as Fco.sMno.sTiO^iSi and 
Euo.4Sro.6S^ In these conventional spin glasses, x" in- 
creases with increasing / both above Tg and below Tg. 

We consider the validity of a dynamic scaling law for 
Tx", which is predicted to be described by24 
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and (c) / dependence of x" below Tg and above Tg. Tg — 
3.92 K. H — Oe. The solid lines are guides to the eyes. 



Tx" = LoynicoT), (3) 
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FIG. 8: H-T phase diagram, where the peak temperatures 
of xzFC vs T (•), dx"{T,HJ = 0.m2)/dr vs T (A), 
x"{T,H,f = 0.1 Hz) vs T (a), and the local minimum tem- 
perature of d5/dT vs T (O). For comparison, Tf{H,uj) with 
/ = 0.01 Hz (■) (see also Fig.Caa)) and Tf{H) (♦) (see Eq. 
^ and Fig. 1141 for definition) are shown as a function of H. 
The solid lines are the fits of the data (see text). 



= 3.92 K, T* = 5.27 x 10"*^ sec. In this plot the data 
X"(T, 7J = 0, w) for 4 < T < 10 K and 0.01 Hz < / < 1 
kHz are used. It seems that the data points fit well with 
a scaling function Q{ujt), indicating the validity of the 
dynamic scaling law given by Eq. 



FIG. 7: (a) T dependence of Tx" at various /. H ^ 0. (b) 
Scahng plot of (Tx")n (= {Tx' )/{Tx")r^a. as a function of 
ujT. {Tx")max IS the maximum value of Tx" ■ The relaxation 
time T is described by Eq. ((TJ with x = 10.3, Tg = 3.92 K, 
and r* = 5.27 x 10"'' sec. 



where ^{ujt) is a scaling function of lot and is assume to 
take a maximum at lot — constant. The value y{= (3/x) 
is a critical exponent, where /3 is a critical exponent of the 
order parameter (q) . FigureUfa) shows the T dependence 
of Tx" at various /. The curve of Tx" vs T exhibits a 
peak, which shifts to the high-T side as / increases. The 
peak height of Tx" defined by {Tx")max increases with 
increasing /. The least squares fit of the data for the 
peak height of Tx" vs uj (for 0.01 < / < 1000 Hz) to 
the form of (« loV) yields y = 0.035 ± 0.001. Then the 
value of (3 {— xy) is estimated a.s (i — 0.36 ± 0.03, where 
X = 10.3 ±0.7. This value of (3 is smaller than that 
{(3 = 0.54) for Feo.sMno.sTiOa-S^ Figure Hb) shows the 
scaling plot of {Tx")n = {Tx") / {Tx")max as a function 
of WT, where r is given by Eq. ^ with x = 10.3, Tg 



D. Possibility of AT line determined from d5/dT vs 

T 

According to the AT theory^^ it is predicted that an 
ideal SG system in the presence of H undergoes a SG 
transition at the spin freezing temperature Tf{H) in the 
{H,T) plane. The deviation of xzfc from xfc starts 
to occur for T < Tf{H). Experimentally it is a little 
difhcult to determine exactly the line Tf{H) at which 
^ — Xfc — Xzfc = because of possible distribution 
of Tf {H) arising from the islandic nature of the system. 
Instead, we use another definition for the line Tf{H) at 
which dJ/dT vs T exhibits a local minimum at each H 
(see Fig. I^b)). In Fig. |H1 we show the hne Tf{H) thus 
obtained (the local minimum temperatures of d5/AT vs 
T) in the {H, T) plane. For comparison, we also plot the 
peak temperatures of xzfc{T, H) vs T, x"iT,H,uj) vs 
T at / = 0.1 Hz, and the local minimum temperature of 
dx"{T,H,u})/dT vs T with / = 0.1 Hz as a function of 
H in the (H, T) plane. Note that the data of x"{T, H, uj) 
vs T at / = 0.1 Hz will be given later (see Fig.l^)). These 
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FIG. 10: (a) - (f) T dependence of x' at / = 0.01 and 0.1 
Hz {H = 0), x' at / = 0.1 Hz at xzfc, Xfc and 5 
XFC - Xzfc) at the same H. H = 100, 200, 300, 500, 700, 
and 1000 Oe. 



FIG. 9: T dependence of (a) x' and (b) x" at various H. f 
= 0.1 Hz and h = 0.5 Oe. 



E. Effect of H on the SG transition 



lines are away from the line Tf{H). The least squares fit 
of the data of the line Tf{H) for 100 Oe < < 1 kOc to 

H = Ho[l-Tf{H)/Tgr, (4) 

yields parameters p = 1.52 ± 0.10 and Hq = (3.4 ± 0.4) 
kOe, where Tg = 3.92 ± 0.11 K. In the mean field 
picture^i*^ this line corresponds to the AT line. In fact, 
the value of exponent p is close to that [p — 1.50) for the 
AT hue. 

It has been believed that the mean field picture is 
vahd for Feo.sMno.sTiOs: p = 1.49.3& However, Matts- 
son et al.A have claimed that the SG transition is de- 
stroyed by the application of a finite magnetic field H 
in Feo.sMuo.sTiOs. Their result supports the prediction 
from the scaling picture that the long range SG order at 
equilibrium only occurs for H = and T < Tg. Further 
discussion for the AT line will be done in Sees. IIII El and 
ElFl 



The transition from the paramagnetic phase to the SG 
phase occurs in our system at H = 0. What happens 
to this transition in the presence of H? Does the long- 
range SG phase survive as suggested by the molecular 
field picture or is destroyed as suggested by the scaling 
picture? According to the method developed by Matts- 
son et al.pi here we have determined the dynamic freezing 
temperature Tf(H, lu) in a wide frequency and field inter- 
val. Figure ini shows the T dependence of x'{T,H,uj) and 
x" {T, H,llj) at various _ff for / = llj/2tt = 0.1 Hz, where 
h = 0.5 Oe and H {0 < H < 3 kOe) is applied along 
the c plane perpendicular to the c axis. In Figs. llOf al - 
(f), we show the T dependence of x'{T, H = Q,uj) dX f — 
0.01 and 0.1 Hz, x'{T,H,uj) at / = 0.1 Hz, xfc{T,H), 
XzFc{T,H), and 6 = xfc {T,H) - xzfc{T,H) at var- 
ious iJ = 100 - 1000 Oe. The susceptibility xfc{T,H) 
{= Mpc/H) corresponds to the equilibrium suscepti- 
bility. Since Mpc{T, H) is a non-linear function of H, 
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FIG. 11: T dependence of x'iT,H = 0,u;) and 
AMfc{T,H)/AH. The spin freezing temperature Tf{H,Lo) 
is defined by a temperature at which x'{T, H = 0,ui) is equal 
to dMFc{T,H)/dH. 



XFC (= Mfc/H) is not equal to the differential FC 
susceptibility dMpc /dH . The numerical calculations 
indicates that dMpc /dH is smaller than xfc at the 
same H and T. In Fig. we show the T dependence 
of dMFc{T,H)/dH at various H and x'(T,iJ = 0,cj) 
at various /. The spin freezing temperature Tf(H,uj) 
is defined as a temperature at which dMpc{T, H)/dH 
coincides with x'{T,H = 0,uj). As shown in Fig. 1101 
X'{T, H, at / = 0. 1 Hz deviates from H ^0,lu) 
at the same / above Tf{H,uj). This is the same method 
that has been used by Mattsson et ali^ to obtain Tf{H, ui) 
for Feo.sMno.sTiOa. The advantage of this method is 
that we do not have to measure the T dependence of 
x' and x" at each H such as the data shown in Fig. |51 
In Fig. I12f a) we show a plot of the line Tf{H,uj) in the 
{H, T) plane, where / (0.01 < / < 1000 Hz) is varied as a 
parameter. For comparison, the line Tf{H, w) at / = 0.01 
Hz is also plotted in the H-T diagram of Fig. |S| We find 
that this line is close to the line Tf{H) determined from 
the local minimum temperature of dS/dT vs T. This 
result suggests that the method to determine Tf{H,uj) 
is appropriate. Note that similar data of Tf[H,uj) has 
been reported by Mattsson et al. for Feo.sMno.sTiOaiS 
Figure I12r b) shows the T dependence of t(T, H) which 
is derived from Fig.^a); t{T,H) = 1/u; = 1/(271 f) at 
T^TfiH,u;). 

Assuming a SG phase transition at Tg, we assume the 
following dynamic scaling relation for the relaxation time 
T±{T, H) in the presence of i?,^^ 

T±{T,H) ^ e F±{X) = i?-2^/('^+'')G±(X) (5) 
with X =1 e I iJ-2/(/3+7)^ where (+) denotes T > Tg 
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FIG. 12: (a) Tf{H,u}) contours at various H between / = 
0.01 Hz and 1 kHz, which is determined from Fig. ill! (b) r 
vs T at various H, which is derived from (a). The lines are 
guides to the eyes. 



and (— ) denotes T < Tg, e is the reduced temperature 
defined hy e = T/Tg — 1, Tc ~ tq \ e is the relax- 
ation time for T > Tg at -ff = 0, 7 is the exponent of 
nonlinear dynamic susceptibility X3i -P+ and are the 
scaling functions for T > Tg and F_ and G_ are ones 
for T < Tg. The scaling relation Eq. (jSJ suggests that 
T diverges like t(T = Tg,H) « iJ-2^/(/3+7) as H ^ 0. 
The characteristic field denoted hy H* «| e |('^+t)/2 is 
a crossover line between a weak-field and a strong-field 
region for T > Tg. 

In the scaling picture, the asymptotic form for t_ be- 
low Tg is obtained as follows. In Eq. Q we assume that I 
diverges like / w iy~2/(d-2e) jj reduces to zero?^ Then 
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T-_ can be rewritten as 

just below Tg, where 9 is the energy exponent. The 
equivalence of Eqs. (|SJ| and (|SJ| leads to the relation 
1/ = (/3 + - 261). Since /3 + 7 = v{d - 29) = 

2 — a — 21^6* = 2/? + 7 — 2^6* , we obtain a scaling relation 
(/3 = 2i/6'). Then in Eq. © has the asymptotic 

form of F_{X) « exp(X'^''), since ^0/(^5^3) w 1. 

Our result of r(T, i?) shown in Fig. I12r b') is ana- 
lyzed using the dynamic scaling relation. The relaxation 
time t(T — Tg,H) is predicted to be proportional to 
^-2x/(/3+7) as i? 0. In fact, the least squares fit of 
the data of r vs at Tg = 3.92 K for 0.1 < i7 < 1 kOe 
yields the parameters '2.x/{l3 + 7) = 5.30 ± 0.13. Since 
X — 10.3 ± 0.7 and /3 — 0.36 ± 0.03, the exponent 7 is 
estimated as 7 = 3.5 ± 0.4. Using the scaling relations 
(a + 2/3 + 7 = 2 and 2 - a — dv), the exponents a and v 
are calculated as a = — 2.2±0.5 and v — 1.4±0.2, where 
d = 3. The exponent z is given byz = 6.6±1.2. These re- 
sults agree with those predicted from the MC simulation 
by Ogielski for the 3D ±J Ising model (a = — 1.9 ± 0.3, 
/3 = 0.5, 7 = 2.9 ± 0.3, 1/ = 1.3 ± 0.1, z = 6.0 ± 0.5)3 
The exponent 9 is estimated as 6 = 0.13 ± 0.02 using the 
scaling relation {6 = (3/2v). The validity of this 9 will be 
discussed in Sec. IIVI 

Figure[T3i;a) shows a scahng plot of t±(T, iy)iy2x/(/3+7) 
a function of X, 

r±(T,i?)i?2-/(/3+7) =G.±(X), (7) 

where x = 10.3, Tg = 3.92 K, /3 = 0.36, and 7 = 3.5. 
We find that almost all the data fall well on two scaling 
functions: G+ for T > Tg and G_ for T < Tg. The 
scaling function F± is related to the scaling function G± 
by F±{X) = X^G±(X) with X = 10.3. In Fig.H^b) we 
show the plot of F±{X) as a function of X. The form of 
is given as follows. ln[F_(X)] = aX"^" + C2 with 
ci = 300 ± 18, C2 = -72 ± 2, and ^'^^ = 0.33 ± 0.01. This 
fitting curve is determined from the data of t_ at -ff = 3 
kOe. Since v ^IA± 0.2, V is estimated as 0.24 ± 0.02. 
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FIG. 13: (a) Scaling plot of t±J/^'=/^''+^' vs X at various H, 
where X =| e | 7J~^/('3+7). The values of exponents (5, 7, and 
X are given in the text, (b) Scaling plot of X'-" t-H^='/^'^+'^^ 
vs X at various H . The solid line is a fitting curve (see the 
text for detail). 



F. Possibility of SG phase transition at finite H 

We consider the possibility that t(T, H) may be de- 
scribed by a critical slowing down given by 

T{T,H)^T*[T/Tf{H)~l]-- (8) 

for T > Tf{H). This is based on the assumption that 
the SG transition is not destroyed by the application of 
H (the molecular field picture). The SG phase exists 
below Tf{H) in thermal equilibrium. The least squares 
fits of the data of t{T, H) vs T shown in Fig. [H to the 
above power law form yield x = 10.2 ± 0.3, Tf{H) — 
3.10 ± 0.05 K, r* = 2.3 x 10"'* sec at iJ = 100 Oe, x 
= 13.2 ± 0.4, Tf{H) = 2.35 ± 0.07 K, r* = 1.2 x lO'^ 



sec at H = 200 Oe, x = 13.8 ± 0.5, Tf{H) = 2.08 ± 
0.08 K, T* = 5.2 X 10^2 sec at H = 300 Oe, and x = 
13.6 ± 0.6, Tf{H) = 1.94 ± 0.08 K, t* = 6.0 x 10"^ 
sec at H = 500 Oe, The values of Tf{H) thus obtained 
are plotted as a function of H in Fig. |S| This line does 
not coincide with the line Tf{H) determined form the 
local minimum temperature of d5/dT vs T. These two 
lines start to deviate at iJ « Oe. The exponent x is 
strongly dependent on H. The value of x at iJ = 100 
Oe is almost equal to x at if = 0, but it drastically in- 
creases with increasing H for H > 200 Oe. These results 
suggest that the SG transition at = is destroyed 
by the application of H (at least above 100 Oe). This 
result is consistent with the following theoretical predic- 
tion proposed by Lamarcq et al.^^ For H < Hq (some 
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FIG. 14: Plot of r vs T obtained from the data of Tf{H,u)) 
with H = 100, 200, 300, and 500 Oe. The relaxation time r 
is defined as lot = 1 at T = Tf{H,uj). 



TABLE L Fit parameters for the decay of x" vs t. T = 3.75 
K to Eq. . The definition of the parameters is given in the 
text. 



/ 

(Hz) 

0.05 

0.1 

0.5 

1 



b" ^iiuj) 

(emu/av mol) 

0.159 ± 0.022 0.522 ± 0.005 0.535 

0.158 ± 0.019 0.502 ± 0.004 0.512 

0.185 ± 0.021 0.445 ± 0.002 0.449 

0.193 ± 0.033 0.419 ± 0.002 0.420 



0.188 ± 0.016 
0.152 ± 0.010 
0.111 ± 0.010 
0.095 ± 0.016 



characteristic field) the system behaves qualitatively just 
as in the case H — 0, while significant changes arise for 
H > Hq). Our result is slightly different from the con- 
clusion derived by Mattsson et aliS that the SG phase 
is destroyed at finite H for Feo.sMno.sTiOa. Note that 
there are some recent theoriesSS supporting the result of 
Mattsson et al. 



G. Aging at constant temperature 

We have measured the t dependence of x" at T = 3.25, 
3.50, 3.75, 3.90, and 4.4 K, where iJ = 0. The system was 
quenched from 10 K to T (< Tg) at time (age) zero. Both 
x' and x" were measured simultaneously as a function of 
time t at constant T. Each point consists in the succes- 
sive measurements of five frequencies (/ = 0.05, 0.1, 0.5, 
1, and 5 Hz). Figure shows the t dependence of x" at 
T = 3.75 K for / = 0.05, 0.1, and 1 Hz. The absorption 
x" decreases with increasing T and is well described by 
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FIG. 15: t dependence of x" at T = 3.75 K for / = 0.05, 
0.1, and 1 Hz, where t is the time taken after the sample was 
quenched from 50 to 3.75 K. The solid lines are the fits of the 
data to the power law form. 



a power-law decay 

x"{u;,t)^x'i{Lu)+A"{uj)t-''". 



(9) 



The least squares fit of the data of x"(t^ii) at T = 3.75 
K to Eq. yields fitting parameters listed in Table 
HJ The exponent b" is slightly dependent on /. The / 
dependence of the amplitude A"{uj) for 0.05 < / < 1 Hz 
is described by A"{uj) = A'^uj-^"" with A'^ = 0.142±0.002 
and /i" = 0.225 ± 0.016. The value of /z" is almost the 
same as that of 6", supporting a cj^-scaling relation that 
Eq. @ can be rewritten as 



X"{u:,t) = x'^{u)+A'^{u:t) 



-b" 



(10) 



In Fig. lTCT a) we show a scaling plot of Ax"(w, t) at 3.75 K 
as a function of iut, where Ax"{(^, t) = x"(w, t) — (xo(^))- 
The stationary susceptibility (xo(t^)) fisted in Table ^ 
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and (b) Ax'(t^, t) [= x'('^i ^ (X(')('^))] ^ ^ function of Ljt for 
the data at / = 0.05, 0.1, 0.5, 1, and 5 Hz. The definition of 
stationary susceptibility (xo(<^)) and (xo('^)) is given in the 
text. The solid lines are the fits of the data to the power law 
form. Each curve is vertically shifted by (xo('^)) or {xo('i')) 
(see text for definition). 





FIG. 17: T dependence of exponents b' and b" 



All the data lie well with a universal curve described by 
Ax'(w,t) = ^fiM)"''' with b' = {b') = 0.226 ±0.007 
and A'f^ = 0.426 ± 0.018. The exponent (6') is nearly 
equal to {b"). The susceptibility (Xo('^)) decreases with 
increasing /, which is similar to the case of (xo(i^))- 

Such a power-law decay of x"{^-,t) ^nd x'i^^'t) with 
t is observed only for 3.25 K < T < Tg. In Fig. [13 
we shows the T dependence of b' and b" at / = 0.05 
Hz. Both b' and b" are on the same order at the same T. 
There is a step-like increase in b' and b" with increasing T 
between 3.5 and 3.75 K; 6" = 0.017±0.032 at 3.50 K and 
b" = 0.16 ±0.04 at 3.9 K just below Tg. Similar behavior 
on the T dependence of b" has been reported by CoUa et 
al.'^-^ for a relaxor ferroelectric Pb(Mgi/3Nb2/3)03, which 
is not a SG system. We note that as shown in Fig. Efb) , 
the data of dS/dT vs T for H = 1 Oe exhibits a local 
minimum at 3.5 K, close to a temperature at which b" 
undergoes a drastic increase. 



is slightly different from Xo('-^) ^.nd corresponds to the 
asymptotic /-dependent value so that x"('^jO tends to 
zero in the limit of t — > oo. Almost all the data fall well 
on a single scaling function given by Eq. with b" = 
{b") = 0.255±0.005 and A'^ = (Af,') = 0.239±0.009. The 
value of (&") is a little larger than that of b" determined 
using the data at / = 0.05 Hz. Note that (&") is equal 
to 0.14 ± 0.03 at 19 K for Feo.5Mno.5Ti03 {Tg = 20.7 
K)i2£ The susceptibihty (xo('^)) clearly decreases with 
increasing uj, which is consistent with the prediction by 
Picco et al;^- 

In Fig. ll6r b') we show a scaling plot of Ax'ii^J, t) a,tT — 
3.75 K, where Ax'(u;, t) = x'(w, t) - {x'oi^))^ and (x[,(^)) 
is the stationary part of x'(^jO the limit of i — > oo. 



H. Rejuvenation effect under perturbations 

We also examine the aging effect of x" under perturba- 
tions such as the change of T and H. Figure^|shows the 
t dependence of x"(w,t) under a negative temperature- 
cycle, where / = 0.05 Hz and h = 0.1 Oe. Here our 
system was quenched from 10tor==3.75Kati/ = 0. 
The relaxation of x"{^^ t) was measured as a function of 
t during a period t^i (« 1.5 x 10^ sec). The tempera- 
ture was then changed to T - AT (the first T-shift). The 
relaxation of x" was measured as a function of t for a 
period t^2 (« 1-5 x 10"* sec) at T - AT. The system 
was again heated back to T (the second T-shift). These 
processes were repeated subsequently. In the case of AT 
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FIG. 18: Relaxation of x"('^i during a negative temperature 
cycle. / = 0.05 Hz, T = 3.75 K and T - AT. H = 0. The 
change of T with f is also shown, (a) AT — 0.05 K. (b) and 
(c) AT = 0.20 K. 



= 0.05 K (Fig.CHIa)), just after the first T-shift x" be- 
haves as if the system were quenched to T — AT directly 
from 10 K. After the second T-shift, however, x"(a;, t) co- 
incides with a simple extension of x"i^,t) already aged 
by t^i at T. In the case of AT = 0.20 K (Figs. HWh) 
and (c)), x" undergoes a drastic change to a value higher 



CUqjC OggDlj-FeClg GBIC 



_l I I I 1 1 I I 1 1 I 1 1 I 1 1 I I I I 1 1 I I I I I 1 1 1 1 I 1 1 I 1 1 I I. 



0.61 



D.6 



0.59 L^H = Oe 



0.58 L 



0.57 




f = 0.05 hfe 

T = 3.75 K 
H = 0- lOOe 



(a) 



jn 1 1 u 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 1 1 1 1 1 u n I n 1 1 . 




FIG. 19: Relaxation of during a magnetic field cy- 

cling. / = 0.05 Hz. T = 3.75 K. (a) H ^ and 10 Oe. (b) 
H = and 50 Oe. 



than the value just before the change of T from 3.55 to 
3.75 K, when T is changed from 3.75 to 3.55 K. Then the 
relaxation of %" newly occurs (the rejuvenation effect). 

Figure [TOI shows the t dependence of x"(a;,t) at T = 
3.75 K under a positive i?-cycle, where / = 0.05 Hz and 
h = 0.1 Oe. In these measurements, the system was 
cooled down to 3.75 K at _ff = 0. The relaxation of x" 
was measured as a function of t during a period t^i. A 
magnetic field H (= 10 and 50 Oe) was then applied at 
t = tiui- The relaxation of x" was measured as a function 
of t for a period t^2 at H. The field H was again reduced 
to zero. The relaxation of x" was measured as a function 
of t at H = 0. These processes were repeated. In the 
case of = 10 Oe (Fig. Efa)), when the field is off 
from 10 to Oe, x" rises to a value of x" to which x" 
is supposed to decay during the period t^2 at H — 0. 
In the case oi H = 50 Oe (Fig. Iiyr b)). x" undergoes a 
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drastic change to a value higher than the value just before 
the change of H from to 50 Oe, when the field is off 
from 50 to Oe. Then the relaxation of x" newly occurs 
(the rejuvenation effect). This effect is similar to those 
reported in Feo.5Mno.5TiO3.a2iS and CdCri. 7^0.384.^° 

IV. DISCUSSION AND CONCLUSION 

In this paper we study the nature of the slow dynamics 
of the short-range Ising SG, Cuo.5Coo.5Cl2-FeCl3 GBIC. 
Using the concepts of static and dynamic scaling laws, we 
determine the critical exponents of this compound. Our 
results are as follows: a = -2.2 ± 0.5, (3 = 0.36 ± 0.03, 
7 = 3.5±0.4, V = 1.4±0.2, ?7 = -0.5±0.2, z = 6.6±1.2, 
</> = 3.9±0.4, 9 = 0.13±0.02, and V = 0.24±0.02, where 
2 — 7y = l/v. These critical exponents are compared 
with those reported for the Ising SG, Feo.5Mno.5Ti03 
(a = -3, /3 = 0.5 ± 0.2, 7 = 4.0 ± 0.4, v = 1.7, 
z ^ 6.2, (/) = 4.5, e « 0.2, and V = 0.3 - 0.7)i9i25,30,33 
and Heisenberg-like SG, CdCr1.7Ino.3S4 {a ~ —1.9, (3 — 
0.75 ± 0.1, 7 = 2.3 ± 0.4, 1/ = 1.26 ± 0.2, = 3.1 ± 0.5, 
zi^ — 7 , = l.l)ffii2i. Our critical exponents are in ex- 
cellent agreement with those of amorphous metallic SG, 
(Feo.i5Nio.85)75Pi6B6Al3 {a = -2.2, (3 = 0.38, 7 = 3.4, 
ly = 1.39, z = 5.9)1^^1^^ Note that we use the scaling re- 
lation given hy (3 = 20v to determine of our system. 
This relation is derived in the present work using the idea 
of dynamic scaling law for the T and H dependence of 
T_. This relation is valid for Feo.5Mno.5Ti03: the value 
of 9 {— l3/2v — 0.15) is close to the experimental value 
[9 K, 0.2). We note that a relation (1 = 2v9) is derived 
by Fisher and Huse?^ in the case of c? ^ . Here df 
is the lowest dimension where 6* > and d'^ is close to 
2.5i2i If (3 is equal to the mean-field exponent (/3 = 1), 
then the scaling relation (/3 = 2v9) coincides with the 
relation derived by Fisher and HusejS^ 

Our value of 6* is a little smaller than the theoret- 
ical values: 9 = 0.19±0.01 by Bray and Moore 
9 = 0.192 ±0.001 by Hartmann,^^ and 9 = 0.20 ±0.03 by 
Komori et al.^^- from MC simulations on the relaxation of 
energy of the 3D Gaussian SG model with nearest neigh- 
bor interactions. Our value of '0 is a little larger than that 
of 0, satisfying the imposed inequality 9 < ip < {d — 1)^ 
This condition is also satisfied in the other SG systems: 
Feo.5Mno.5Ti03 (V- = 0.3 - 0.7)^ and CdCr1.7Ino.3S4 
{ijj = l.l)i2& The value of ip is strongly dependent on the 
spin symmetry (n) of the systems such as Ising {n = 1) 
and Heisenberg (n = 3): it decreases as the spin symme- 
try n decreases. This is also true for the critical exponent 
/3: P decreases with decreasing nM- The small values of 

and (3 indicates that our system magnetically behaves 
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like an ideal 3D short-range Ising SG. 

Finally we discuss the T dependence of b". The expo- 
nent b" is dependent on T. According to Komori et al.^ 
the ratio Ax"(w,t)/x'(^) is described by 

Ax"{cu,t)/x"{io)»{ojt)-''"^^\ (11) 
where 6"(T) = (d - 9)/z{T) and the exponent l/z{T) 
is linearly dependent on T well below Tg-. l/z{T) = 
T/{zTg). This wi-scaling of Ax" (uj , t) / x" (uj) is derived 
on the assumption that the growth law of the spin- 
glass correlation length ^(t) is approximated by the form 
^{t) lo{t/toy^^^^\ where Iq and to are microscopic 
length and time scales. This growth law is different from 
that proposed in the scaling picture due to Fisher and 
Huse.^ Using the values of 9 and b"{T), the exponent 
l/z{T) is calculated as 0.006 at 3.50 K and 0.056 at 3.90 
K. Such an increase of 1 / z{T) with T, qualitatively agrees 
with the prediction by Komori et al42 If the expression 
l/z(T) = T/{zTg) is valid at T = 3.90 K just below Tg, 
the value of z can be estimated as 18, which is much 
larger than our value z (= 6.6). Experimentally b"{T) 
is nearly equal to zero below 3 K (« 0.75Tg). At the 
present we give no satisfactory explanation for the cause 
of the drastic change in b" around 3.5 K in our system. 

In conclusion, we have shown that Cuo.5Coo.5Cl2- 
FeCl3 GBIC magnetically behaves like an ideal 3D short- 
range Ising SG. This compound undergoes a SG phase 
transition at Tg — 3.92 ± 0.11 K. The dynamic scaling 
analysis suggests that this SG transition is destroyed by 
the application of H (at least above 100 Oe). The scaling 
behavior of the relaxation time is well described by the 
scaling picture with the energy exponent 9 = 0.13 ± 0.02 
and the barrier exponent -0 = 0.24 ± 0.02. The aging 
obeys the wi-scaling. The rejuvenation effect is observed 
under sufficiently large (temperature and magnetic-field) 
perturbations. 

In spite of such a similarity in dynamic behaviors 
between Cuo.5Coo.5Cl2-FeCl3 GBIC and 3D Ising spin 
glass, the frequency dependence of x" below Tg are 
very different. The absorption x" in Cuo.5Coo.5Cl2-FeCl3 
GBIC decreases with increasing frequency, whereas x" in 
3D Ising spin glass increases with increasing frequency. 
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